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I . INTRODUCTION 



One of the most important issues in the world today is 
the soaring prices of energy resources, this has caused con- 
siderable anxiety in the energy consuming nations concerning 
their ability to maintain economic growth and financial 
stability. 

For example, in his second economic report, President 
Nixon cautioned that: 

"The new problems we face are of such enormity 
that there may be a temptation to delay further 
progress toward trade and monetary reform. The 
drastic increase in oil prices will have a signif- 
icant impact on both the domestic economies of all 
nations and international economic relationships . [1 ] 

According to United Nations Statistical Office, the total 
energy consumption in the world increased by 4% in 1971. 

The world consumption of electricity increased by 6% over 
previous year; world production of crude petroleum increased 
by 5.5%; the world exports of natural gas increased by 36% 
over 1970. The nuclear power accounted for 2% of total en- 
ergy produced, which was an increase of 38% over 1970. On 
the other side the use of solid fuels declined by 9% of 
total world energy consumption. [2] 

Government revenues in oil producing nations have in- 
creased considerably due to the recent increases in the price 
of crude oil. In fact the posted price of Persian Gulf light 
crude oil was increased from U.S. $1.80 per barrel on June 1, 
1973 to U.S. $11.50 per barrel on January 1, 1974 by OPEC 
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(organization of Petroleum Exporting Countries), consequently 
producing nations are accumulating more reserves and their 
surplus could rise from $10 billion in 1973 to $60 billion 
in 1974. [3] 

On the other side, consuming nations, due to the higher 

oil prices would have to have current' account deficits. In 

their conference in Washington the consuming nations agreed 

that: "...the present situation, if continued could 

lead to a serious deterioration in income and 
employment, intensify inflationary pressures, 
and endanger the welfare of nations. And a 
serious setback to the prospect for economic 
development of developing countries will re- 
sult from their additional energy costs at 
current oil prices." [4] 

The production and consumption of energy has significant im- 
pact on the economy and financial position of a country. 

Hence in the last few years, many studies have been conducted 
to evaluate the impacts of policy options such as price reg- 
ulations, import regulations, environmental controls, etc. 
and several of these studies have used advanced analytical 
techniques . 

In order to be able to trace the impact of policy changes 
and to be able to develop an analytical model for the con- 
junctive management of the energy sector of economy, it is 
essential to have an understanding of some of the complex 
interrelationships in the broad energy, demand and pricing 
situation . 

This is the primary objective of this study. By conduct- 
ing a conprehensive analysis of the world's energy consump- 
tion, energy-GNP relationship, and the existing literature, 
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an analytical model has been developed. The usefulness of 
this study for future efforts, in analyzing the consequences 
of alternative energy policy options, is covered in the con- 
clusion. 
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II. AN EMPIRICAL STUDY OF THE ENERGY-GNP RELATIONSHIP 



A. INTRODUCTION 

"Economic history attests to the crucial role 
played by the consumption of energy in advancing 
the material well-being of mankind. The industrial 
revolution and the growth of industry in the nine- 
teenth century are almost synonymous with the sig- 
nificant contribution of coal to the development 
of the iron and steel industry, to railways, and 
to factory mechanization. In the twentieth century, 
electrification and motorized transport served to 
modify the process and many ways helped economic 
progress. 'High-energy civilization' has by now 
become a commonplace characterization of economi- 
cally-advanced 20th Century society. The growth 
of population and the levels of prosperity achieved 
have depended upon, and resulted in (among other 
critical factors) the consumption of energy in 
large amounts. To the extent it materializes, con- 
tinued demographic and economic growth will un- 
doubtedly involve the steadily rising use of ener- 
gy." [5] 

It therefore becomes important to give particular attention 
to this relationship. At a minimum, the topic calls for an 
interlinking of physical energy data and of the nations' in- 
comes and products (GNP) accounts. 

Such a study could yield a number of important results: 

A better understanding of the behavior of energy 
consumption relative to gross national product. 

It would attempt to answer the question: On the aver- 
age what levels of energy consumption are required 
for given levels of GNP. [6] 
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B. DEFINITION OF THE MODELS 



We now turn to the specification of a basic model relat- 
ing GNP and energy consumption. 

• In the first model the assumption is made, that a linear 
relationship exists between per capita energy consumption 
and GNP per capita, then: 

(1) E = a+BG+e 

where; E = energy per capita, (Dependent variable) 

G = GNP per capita, (Independent variable) 
a = constant 

B = the coefficient representing the propensity to 
consume energy, (in physical units) 
e = stochastic disturbance. 

Next the assumption is made, that an exponential rela- 
tionship or logarithmic relation exists between these vari- 
ables, then: 

(2) E = yG 6 I , or 

(3) In E = lny + 51nG + lne 
where; y = lny = constant. 

6 = coefficient representing the elasticity of energy 
consumption with respect to GNP 
e = lne = stochastic disturbance 
The models represented by equations (1) and (3) were applied 
for a number of samples of data corresponding to OECD Coun- 
tries, (Organization for Economic Cooperation and Develop- 
ment). The results of the regression analysis are presented 
in section D, and the summary of principal statistical re- 
sults appear in Table I. 
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C. DESCRIPTION OF DATA 



The data used here are collected from United Nations 
Statistical yearbooks and Statistics of Energy for member 
countries of OECD listed in Appendix C. The selection of 
data was made because of the special nature of the problem 
which requires detailed statistics of GNP in constant prices, 
with corresponding values for population and energy consump- 
tion. The only available data with the appropriate charac- 
teristics were those of OECD countries and only for selected 
years, the specific data which is used is described in fol- 
lowing lines. 

1 . Population 

Data are taken from OECD Main Economic Indicators, 
and is measured in million, for years 1960, 1965, 1966, 1967, 
1968, 1969, 1970 and 1971, total of 174 observations (22 
countries multiplied by 8 years minus two observations for 
years 1966, 1967 of Australia, not being used, because of 
lack of corresponding values for energy consumption per capi- 
ta during these two years. ) 

2 . Energy Consumption per Capita 

Data are taken from United Nations Statistical year- 
books; 1960, 1965, 1966, 1967, 1968, 1969, 1970 and 1971. 

It is measured in kilograms coal equivalent per capita. 

There are 174 observations. 

3 . Gross National Product 

Data are taken from OECD Main Economic Indicators. 

GNP is measured in million U.S. dollars, at 1963 prices and 
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1963 exchange rates. The GNP figures are next divided by 
corresponding values of population to determine the GNP/ 
capita, for the years 1960, 1966, 1967, 1968, 1969, 1970 and 
1971 with a total of 174 observations. 

4. By_ Sources of Energy 

All data are taken from "Statistics of OECD countries 
for Energy, 1954-1968, Paris 1970" for the years 1966, 1967 
and 1968, and divided by corresponding population figures 
(Appendix C). 

a. Crude Petroleum 

This is measured in thousands metric tons with a 
total of 60 observations (20 countries for the three years, 
1966, 1967 and 1968). 

b. Motor Gasoline 

This is measured in million metric tons with a 
total of 63 observations (21 countries for the three years, 
1966, 1967 and 1968). 

c. Aviation Fuel 

This is determined by sum of the aviation gaso- 
line and jet fuel. It is measured in million metric tons 
with a total of 63 observations (21 countries for the three 
years, 1966, 1967 and 1968). 

d. Hard Coal 

This is measured in million metric tons with a 
total of 63 observations (21 countries for the three years, 
1966, 1967 and 1968). 



14 



e. Kerosine 



This is measured in million metric tons with a 
total of 63 observations (21 countries for the three years, 
1966, 1967 and 1968). 

f. Gas Oil 

This is measured in million metric tons with a 
total -of 63 observations (21 countries for the three years, 
1966, 1967 and 1968). 

g. Natural Gas 

This is measured in billion cubic meters at 4200 
KCal, with a total of 33 observations (11 countries for the 
three years, 1966, 1967 and 1968). 

D. STATISTICAL RESULTS 

For the application of the models a general-purpose sta- 
tistical package, based on least-square method (SNAP/IEDA) 
has been used. The results of the analysis are as follows 
(with standard errors in parenthesis), 
where; R 2 = Square of correlation coefficient. 

R 2 = Adjusted R 2 . 

F = The ratio of the variance of the residuals of 
the dependent variable before regression and 
the variance of the residuals of that variable 
after the regression. 

NOBS = Number of observations. 

M = The mean of variables. 
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K = Elasticity of consumption. When the linear mod- 



JE MG 
3G ME' 



el is used K = -£ 7 ? 77 ^, where is the regres- 



sion coefficient, MG and ME are the means of 
the independent variable and dependent variable 
respectively. When the constant elasticity 
model is used, K equals the coefficient 6. 
Energy/Capita and GNP/Capita 

a. Linear Model (Fig. 1) 

E = -584.905 + 2.625G 

(.093) 

R 2 = .823, F = 797.1, NOBS = 174, 

ME = 3722.8621, MG = 1641.6658, 

K - i'l - 2 - 625 * iftBS? - 2 - 625 * - 441 

K = 1.158, R 2 = .822 

b. Constant Elasticity Model (Fig. 2) 

InE = -0.902 + 1.225 InG 

( .031) 

R 2 = .903, F = 1603.4, NOBS = 174, 

K = 1.225, R 2 = .902 

Crude Petroleum/Capita and GNP/Capita 

a. Linear Model (Fig. 3) 

Ep = 269.059 + 0.709G 

(0.269) 

R 2 = .107, F = 6.9, NOBS = 60, 

MEp = 1401.8915, MG = 1597.0630, 

1597.0630 

K - -709 x 14Q1 g915 

K = .709 x 1.14 = .81, R 2 = .092 
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b. Constant Elasticity Model 



(Fig. 4) 



InEp 


= -0.334 + 0.998 InG 
( .126) 


R 2 = 


.52, F = 62.8, NOBS = 60 



K = .998, R z = .488 

3 . Motor Gasoline/Capita and GNP/Capita 

a. Linear Model (Fig. 5) 



Eg = 


-172.334 + 0.258G 
(.021) 


II 

CM 

P5 


.705, F = 145.9, NOBS = 63, 


MEg = 


= 244.5014, MG = 1616.1303, 


K = . 


_ .. 1616.1303 
258 244.5014 


K = . 


.258 x 6.610 = 1.705, R 2 = .700 



b. Constant Elasticity Model (Fig. 6) 

InEg = -4.421 + 1.321 InG 

( .061) 

R 2 = .885, F = 470.6, NOBS = 63, 

K = 1.321, R 2 = .883 
4 . Aviation Fuel/Capita and GNP/Capita 

a. Linear Model (Fig. 7) 



Ea = 


-23.644 + 0.045G 
( .006) 


II 

CM 

cn 


.451, F = 50.1, NOBS = 63, 


MEa = 


= 48.4012, MG = 1616.1303, 


K = 


„ 1616.1303 
■ 045 * 48.4012 


K = 


.045 x 33.390 = 1.503, K 2 - .442 
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b. Constant Elasticity Model 



(Fig. 8) 



InEa = -4.304 + 1.061 InG 

( . 166 ) 

R 2 = .402, F = 40.9, NOBS = 63, 
K = 1.061, R 2 = .392 
5. Hard Coal/Capita and GNP/Capita 
a. Linear Model 



b. 



Eh = 
R 2 = 
MEh 

K = 

K = 



69.885 + 0.470G 
( .123) 
.192, F = 14.5, 



828.7752, MG = 



470 x 



1616.1303 

838.7752 



470 x 1.95 = .92 



NOBS = 63, 
1616.1303, 



R 2 = .179 



Constant Elasticity Model 



InEh = -.917 + . 957G 
( .267) 

R 2 = .174, F = 12.9, NOBS = 63, 
K = .957, R 2 = .160 



(Fig. 9) 



(Fig. 10) 



6 . Kerosine/Capita and GNP/Capita 



a. 



b. 



Linear Model 



(Fig. 11) 



Ek = 
R 2 = 
MEk 
K = 

K = 



6.973 + . 014G 
( .005) 

.121, F = 8.4, NOBS = 63, 



= 29.6929, MG = 1616.1303, 



.014 x 



1616.1303 

29.6939 



.014 x 54.43 = .762/ R 2 = .107 



Constant Elasticity Model (Fig. 12) 

R 2 = .06, R 2 = .045 
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7 . Gas Oil/Capita and GNP/Capita 



a. Linear Model (Fig. 13) 



Eo = 


= 229.796 + . 452G 
( .065) 


R 2 = 


= .443, F = 48.6, NOBS = 63, 


MEo 


= 960.6031, MG = 1616.1303 


K = 


.. 1616.1303 
* 960.6031 


K = 


.452 x 1.682 = .760, R 2 = .434 



b. Constant Elasticity Model (Fig. 14) 

InEo = -1.462 + 1.121 InG 





(.071) 


R 2 = 


= .805, F = 251.7, NOBS = 63, 


K = 


1.121, R 2 = .802 


8 . Natural 


Gas/Capita and GNP/Capita 



a. Linear Model (Fig. 15) 



En = 


= -2013.059 + 1.708G 
(.184) 


R 2 = 


= .736, F = 86.4, NOBS = 33, 


MEn 


= 992.0987, MG = 1759.8145, 


K = 


1 70 o v 1759.8145 
1,708 • 992.0987 


K = 


1.708 x 1.77 = 3.03, R 2 = .727 



b. Constant Elasticity Model (Fig. 16) 

InEn = -27.099 + 4.392 InG 





( .609) 


R 2 = 


= .626, F = 52.0, NOBS = 33, 


K = 


4.392, R 2 = .614 



The results of the empirical analysis are summarized in 
Table I. 
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SUMMARY OF THE RESULTS 
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E. INTERPRETATION OF THE RESULTS 



By referring to the summary of results presented by Table 
I, we can reach the following conclusions: 

1) The logarithmic model was a better fit to the data 
on energy per capita, crude petroleum, gasoline, and gas oil. 
While the simple linear model was a better fit to data on 
aviation fuel, hard coal, kerosine and natural gas. 

2) There is a higher correlation between GNP and aggre- 
gate energy consumption than between GNP and different ener- 
gy sources. This could be explained by the fact that each 
member country of the OECD consumes different sources of 
energy as the primary mixture for energy production. There 
are a number of economic and technological factors affecting 
the mix of the primary energy sources within a country. It 
is therefore necessary to pay particular attention to the 
special factors underlying the historic interfuel composition 
of total energy consumption and to the policy tools avail- 
able to alter the composition. 
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III. THE ANALYTICAL STUDY OF ENERGY PRODUCTION AND MANAGEMENT 



A. INTRODUCTION 

There are considerable efforts taking place in the private 
and public sectors of the economics profession in the recent 
years, to satisfy the need for a quantitative economic theory 
of energy conservation and management. Such efforts provide 
a point of departure for the development of analytical mod- 
els which may be usefule in projecting future energy require- 
ments and the short-term, and the long-term adequacy of 
supplies to meet them. 

The energy situation comprises a wide range of interre- 
lated socio-economic and institutional factors. These in- 
clude: 

- Balance of payments. 

- National security and stockpiling. 

Environmental costs and user costs. 

- Protection of energy related industries. 

Resource production (exploration, development, ex- 
traction, processing, refining) and relative facil- 
ities . 

- Transportation (means, modes and facilities). 

- Marketing and Distribution systems. 

Substitutabilities among different energy sources. 
Development of new sources of energy. 



22 



- Problem involving risks and uncertainties in terms 
of domestic and foreign investments in exploration 
as well as for import supplies, prices, and uncer- 
tainties related to technologies. 

Interdependencies among the stages of activity and 

4 

the policy alternatives. 

Conjunctive management at national, regional, and 
state level versus an individual firm or a single 
industry . 

The purpose of this chapter is to present a summary of 
the literature concerning the production and management of 
the natural resources. 

B. RESULTS OF THE ECONOMIC THEORY OF EXHAUSTIBLE RESOURCES 

The theory of exhaustion considers the losses firms may 
incure if they do not consider the resource availability in 
future generations. 

Results indicate that different structure of industry may 
produce different output patterns, for example private owner- 
ship may lead to decisions different from those of central- 
ized management. In addition, the time-path of exploitation, 
when the costs are affected by cumulative production is dif- 
ferent from conditions where the increasing costs due to 
cumulative output are ignored. 

The general format of the models presented under Exhaus- 
tion theory consist of maximization of present value of 
future profits subject to limited supply of resources. In 
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following section simple examples of these models have been 
presented. 

C. SIMPLE CASES IN THEORY OF EXHAUSTION 
1. Case CD 

Considering the case where a private distributor of 
fuel is hoarding a fixed quantity of it for a variable length 
of time without occurring additional costs, and the selling 
price of fuel increases by a known function P = p(t); as 
time goes on. What would be the optimum time for him to sell 
stored supply knowing that he can sell all quantity Q avail- 
able? [7] 

It is assumed that interest is added continuously at 
a constant rate of r%. 

If (y) denotes the present value of the fuel sold 
after (t) years, then 

(1) y = Q.p(t)*e" rt 

The condition for optimum time are given by: 




Knowing that the log y is a "monotonic transformation" of y. 
(This means that whenever y is increasing log y is also in- 
creasing, and whenever y is falling its logarithm is also 
falling.) Therefore, the point corresponding to the maximum 
of y is also maximum of log y. Then 

(2) log y = log P(t) + log Q - rt, or 

(3) log y = F(t) + log Q - rt 
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where F(t) = log P(t). Then: 



dt (log y) 



1 d£ . 
y dt 



F ( t )-r 



dt 



.(logy) =^:(i If) = i |^ - ir(ff , 2 = f" 



dt v y dt 



y dt 



(t) 



The condition for optimum selling time are thus 
F (t) <0 and F (t) = r 
From Eq. (3) the first condition is that 

f 

F’(t) = ^(log P(t)) = = r 

Where ^ determines the rate of growth^ - of selling price 

at period (t), is equal to the constant rate of market inter- 
est r. By the second condition 



F (t) = dtf^ 1 ^ P(t)) = 7TT{ 



d_ P(t) 



dt l P(t) 



} < o 



which implies the rate of growth for selling price must be 

decreasing at the optimum time. 

Result: If the rate of increase for selling price 

is constant and equal to the market rate of interest over the 

time; then owner would be indifferent whether he sells his 

rt 

reserves at price Pq now or at a price P=PqC after time t. 
This could be shown from Eq. (1). The present value would 
be: y = QPe -rt 

pt 

where P=Pq 6 implies 

(4) y = QP 0 e rt . e “ rt = QP q 



Appendix A, Section 1. 
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2 . Case (2) 



Now considering the case under free competition when 

the rate of increase for selling price is equal to the mar- 

rt 

ket rate of interest. The formula P=P^e fixes the relative 
prices at different times and the owner would be indifferent 
about a unit of fuel being sold now or in the future. (P 
being interpreted as the net price received after paying the 
costs of extraction and placing upon the market.) [8] And Q 
being the total amount of supply available, and the quantity 
q distributed at any given time is a continuous function of 
price and time, if t=T being the time of final exhaustion, 
it could be written: 

(5) q = q(P,t) when t=0, P=Pq; then Q = Qq* 

and when t=T, P=PQe r ^ and q=0. Then 

T T 

(6) Q = / q(P,t)dt = / q n (P n e rt ,t)dt 

0 0 

The question will be the determination of t=T the 
time of final exhaustion to maximize the present value of 
his revenue. Since at t=T, q=0 it implies: 

(7) q = q(P,t) = q(P 0 e rT ,T) = 0 
and consequently 

T T 

y = / q(P, t ) •P*e _rt dt = / q(P , t ) • Pe rt • e _rt dt 
0 0 

T T 

= / q(P, t ) *P dt = P n / q(p, t )dt 
0 0 

From Eq. (6) follows y=PgQ where P Q and Q are both constant. 
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Result: In the case of free competition, when also 

the rate of price increase is constant and equal to market 
rate of interest, the owner is indifferent to q quantity 
supplied at any time t from Q his total reserve and the T 
final time of exhaustion. 

D. A GENERAL CASE IN THEORY OF EXHAUSTION 

Considering the case under monopoly when a certain com- 
modity is available in fixed supply Q with the rate of pro- 
duction q=q(t) over a time period from t^ to [9, 10]. 
Production then is subject to the constraint. 
t 2 

(8) / q(t)dt < Q 
t o 

The firm would plan to maximize the present value of its 
profits 

t2 -rt 

(9) y = / II[q(t),t]*e dt 

to 

where y is net present value; n total profits; q(t) output; 
r the continuous interest rate. 

The question involves the selection of the optimal times 
for starting and completing output and the optimal output 
pattern. The solution implies resolving the Lagrangean Equa- 
tion . 

t 2 _ rt t2 

(10) L = / II [q( t ) , t ] • e r at - X f q(t)dt-Q . 

to to 

where X is the Lagrangean multiplier; and Q total availabil- 
ity. 
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The optimal production pattern is determined by the 

2 

Euler equation of the Calculus of Variation, and the bound- 
ary conditions are the attaining stationary value at t 0 and 
1 2 (Eq. 10). The Euler equation requires that 

9L _ _9 , 8Lv _ j-. 

9q at' 9q ; 

And because does not appear in Eq. (10), it therefore 

does not effect profit, so the Euler equation reduces to 
9L 



9q 



= 0 . [ 11 ] 



/ \ 3L -rt /* 31T 9 q j , ~ 

(11> 15 = e ' 15 lt dt - Xq = 0 

t 0 

By substitution of dq for -|~[dt ; Eq . (10) could be written as 



( 12 ) 



9L _ -rt Van 
9q 



= e 



l 15 dq ' xq or 



(13) / -l-^dq = Xqe rt or 

t 0 3q 

(14) M = Ae rt 

By introducing the concept of marginal profit, the revenue 

(R) minus Cost (C) would determine the profit. Using letters 

9R 



P for price; AC, average costs; <f> , average profits; MR = 



9C 3 n 

marginal revenue; MC = — , marginal cost; and M = — , 

ginal profit. Then Eq. (14) can be written as 

(15) 



9q ’ 
mar- 



Mn(t) = |^ = MR - MC = Ae rt 



The conventional second order equations 

9 2 R „ 9 2 C , 

rr < — r also apply. 



9 2 n 

aq 7 



< 0 and 



9q" 



3q‘ 



Appendix B. 
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In fact II = R - C = P(q)*q - C(q), 



(16) 



Mn=||=||-||= P(q) + q-P(q) - C(q) and, 

j\2tt a 2 P £ 2 p 1 • II II 

4 = 4^-^ = P(q) + P(q) + q-P(q) - C(q) 



9q z 9q z 9q' 
where in general P(q) being the demand curve implies its 

I II 

slope P(q) < 0 and P(q) 2l 0. Considering increasing cost 

II II 

industry, it then implies C(q) > P(q) > 0;’ consequently mak 
ing Eq. (16) negative or 

9 2 n „ „ 9 2 R 9 2 C „ „ 9 2 R „ 9 2 C 

w °' or V "Sq 1 °' or 3^ 3^ • 

Under free competition, P = MR so Eq. (15) becomes 
(17) P - MC = Xe rt 

and with constant cost, this reduces further to 



(18) 



P - AC = Xe 



rt 



Result: From Eq. (15) it derives: 

(19) n = Xq(t)e rt 
dividing it by q(t) results: 

(20) <Kt) = Xe rt 

since Eq. (15) also hold at time t 2 , the optimum production 
path is where marginal profit equals average profit. Since 



II(t) = 4>(t)*q(t), then Mll(t) = 4>(t) + (^)q(t). 

The equality MII(t) = <J> ( t ) implies q(-|~) = 0 which occurs if 
q or is zero; the firm either produces nothing at t 2 or 
produces the output that maximizes the average profits at t 
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E. RESULTS OF ECONOMIC MODELS OF PRODUCTION IN RESOURCE 

INDUSTRIES 

The theory of exhaustion has been developed without ex- 
plicitly considering the role of investment associated with 
resources, the interaction among production rates, invest- 
ment rates, technological changes and exploration expendi- 
tures . 

In recent studies efforts have taken place to develop 
economic models to include these factors. According to these 
studies, the stages of activity necessary for the production 
of energy commodities involve numerous interrelated aspects. 
For example, the optimal extraction rate of a given resource 
and associated investments in capital stock are related to 
exploration and importation policies of the resources as 
well as the substitutability among the resources. [12] 

The supply of particular resource should be obtained 
from domestic production and imports until marginal import 
costs equal the marginal costs of extracting domestic stocks, 
and be allocated between current use and stockpiles for fu- 
ture consumption, until the marginal contribution of resource 
to the current use equals its marginal contribution to future 
use. [13] 

By considering the user costs defined by A. Scott as "the 
present value of future profits foregone by a decision to 
produce a unit of output today." The use of imports in cur- 
rent production would reduce the rates of domestic extraction 
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and consequently the user cost, or imports could be used to 
increase the stockpiles to reduce the future scarcity of the 
resource. [14] 

The above brings to attention the costs which accompany 
the import restrictions and/or stock controls. 

The equality of the user costs specially among substitut- 
able resources are of special interest in the conjunctive 
management of energy-related resources. 

An important question arises in terms of the impact of 
subsidies, depletion allowances granted to one resource, on 
the optimal rates of exploration expenditures, importation 
and production of other resources. [15] 

Of particular attention are the social attitudes towards 
the generation of pollutants, and the inclusion of pollution 
abatement costs in the processing costs and production costs, 
which may be imposed in terms of social costs. 

Considering the problems associated with petroleum pro- 
duction and reservoir management, the results suggest that 
the amount of pressure which exists in the reservoir, and 
the time-path of pressure depletion, are two influential fac- 
tors on the production rate and the amount of oil recover- 
able [16] 

It should be noted that in many cases pressure may be 
maintained or increased artificially by incurring additional 
expenditures, which implies that the recoverable stock, as 
well as production rates, generally will depend strongly on 
investment. This dependence is not merely on a cumulative 
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investment, but also on the time-path of investment. Thus 
the impact of rapid production rates, in the face that the 
recoverable stock is rate dependent, may be offset to some 
extent by investments. [17] 
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IV. DEVELOPMENT OF AN ANALYTICAL MODEL 



A. INTRODUCTION 

Because of increasing concern over the management of 
"Energy Sector" of economy as whole, a comprehensive analyti- 
cal model is required, to be able to translate and quantita- 
tively evaluate the impact of energy related procedures. 

In the three preceding chapters a brief review of litera- 
ture, past and current works was performed, to give an under- 
standing of some of the complex interrelationships in this 
sector. 

In specifying the scope and coverages of this model, a 
number of requirements, could be identified, to be covered 
by the comprehensive model. These requirements could be 
listed as follows: [ 18 ] 

The model should cover all energy forms. 

- The model should consider supply and demand inter- 
actions or the mechanism for balancing supply and 
demand . 

The model should take into account interfuel substi- 
tution . 

The model must take into consideration the impact of 
price and other variables on. the demand and supply of 
the specific fuels. 

The model must consider the impact of technological 
changes on both the supply and demand of fuels and 
energy . 
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The model must be disaggregated at the regional level. 
- The model must be able to use the available data for 
future analyses. 

Note that the environmental consequences have been left out 
in this study. 

The economic model prepared by R. Cummings and David 
Whipple [19] , is an initial attempt at a comprehensive model 
of the discovery, extraction, importation, and stockpiling of 
primary energy resources, integrating the production of final 
energy commodities, it is believed its further expansion and 
refinements, will make it suitable for use in constructing 
an analytical model, presented in next section, feasible to 
the conjunctive management problems of energy sector, by 
covering the above mentioned requirements. 

B . THE MODEL 

By assuming the existence of "j" resource fields in which 
any of "R" resources are extracted from domestic deposits 
and/or stored in stockpiles. For each field j=l,2,...,J; the 
intertemporal changes in stocks of domestic deposits of the 
resource r, r=l,2,...,R, are described by the following: 

(1) X t+1 = X t .-W t ., X°. given, for all r=l,...,R; 
rj r j r j rj 

3=1, • • • , J; t=l , 2 , . . .T. 

By (1), the sotck of resource r in field j at the end of 

period t, "X t+ ^" equals stocks at the beginning of the period 
J 

t t 

M X minus extraction of r during t, "W which is a 
rj b ’ rj 

function of : 
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= The excess of "P^" ; the market price of re- 

wt 

source r, over "C the extraction cost of 

rj 

resource r in field j, during period t. [20] 



. = The excess of P 
rj 



r 



over "c^" the importation 



cost of resource r in field j , during period 



t. [21] 



= The quantity of resource r, demanded in period 
t, which is a function of P^; r=l, . . . ,R the 
market prices of all resources during the 



period t . 




It is alos assumed, that periodic extractions of r in field 
j are subject to an upper bound of the form. 



resource and capital stocks in terms of impinging on periodic 
rates of extraction. For example, given a capital stock 



stock of petroleum in petroleum reservoir (and therefore, the 
lower may be the pressure for natural derive), the samller 
one would expect the maximal periodic rate of extraction. 

For a given resource stock, smaller capital stocks imply 
smaller periodic rates of extraction. 

Further, it is assumed that r extracted from j's deposits 
are used for increasing stockpiles of r. 
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The intertemporal changes in stockpiles of the resource 
r, are described by the following: 

oO 



(5) S t+1 = S 1 .+W t .+M^ .-A t . 
v rj rj rj rj rj 



rj given, for all r, j and t. 



By (5), stockpiles of r, in field j at the end of period t 
equals initial stockpiles, "S^"' plus extraction 
"W^ " , plus imports "M^." , minus substract ions from stock- 
piles "A L " for current consumptions. 

J 

The following functional forms are assumed for A^ and 



M*.. 














rj 














(6) 


A^ 

rj 


= <j>,. .(n't . 

<5j rj=l, . 


t 

..,J’ P rj=l,...,J’ 


y rj=l, . . 


• ,J’ 


Q*). 




9<t> 6 


3 < j ) 6 


3(j) 3(f) 








with 


an 


> °» 3p > 


i °> 3Q 1 0 








(7) 


11*. 
r J 


= <t> .(nt . - 

mj v rj=l, . 


t 

. .,J’ P r j=l , . . . , J ’ 


H 

II 

4-> U 
P- 


. ,J’ 


Q l ) 

r' 



t s t 

where; y . = The excess of P , over "C . " the storage cost of 
rj r’ rj & 

resource r in field j, during period t. 
Stockpiles S . are assumed to be restricted by an upper bound 
which reflects the capital stocks and storage facilities. 
With defined as field's J capital stock for storage fa- 

cilities at the beginning of period t, it is defined. 

(8) . < N*'.(K t .), ^ > 0 for all r, j and t. 

r J — r j v r j ' 3K — 



In the resource extraction system as characterized by 
(1)— (8) , two observations are of interest at this juncture. 
First, statement (5) allows for the inclusion of port facil- 
ities for imports to be treated as "resource fields", i.e., 
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* t t 

for some j=j , X . W . * are zero, and resource management 

r J 

* 

issues in terms of "field" J concerns the management of 
stockpiles of r from imports. Second the transition equation 
for stocks given in (1) is in a simplified form for exposi- 
tory purposes, and does not allow for the dependence of 
recoverable stock on (earlier) time-rate of production and 
investment which may be of particular “importance for re- 
sources such as petroleum. [22] 

The system (l)-(8) is assumed to apply to fields which 
are known (for example, "proven" reserves at t=l) at the 
initial time-period. 

To allow however, for exploration expenditures, e , 

i Z 

which may result in new fields, z-l,...,Z, with deposits of 
the resource. 

Let exploration expenditures e^ z result in discovery of 
h (e ) units of the resource r. (h, is essentially a meas- 
ure of the "probably" rewards of exploration expenditures). 

By assuming a concave set of relationships of the form: 

T t 

(9) £ h rz (e rz ) - X rz’ for a11 r ’ z ’ and 



3h 3h - 3h n 

rz > 0, > r2 



3h 



rZ 



3e 



rz 



3e - 3e „ 
rl r2 



3e 



rZ 



By (9), initial exploration expenditures are for h ^ (which 

yields larger "Probable rewards" per unit of expenditure than 

h rz , z=2,...,Z), until a bound X^ is effective, after which 

h 0 ,...,h „ are applicable. The bounds X are included to 
r ^ rz rz 

allow use of a continuous model as well as to allow for 
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establishment of new fields, (with reserves X ) as a result 

rz' 

of exploration expenditures. It is also assumed e£ z to have 
the form. 

(10) e t = $ (n , p , n , Q„), 

v ' rz e rz rz rz r 



with 



3 <p 



3ii - 



- o, 



3<J> 



3 P - 



: 0 , 



3<f) 



an - 



> 0, and 



3P 



3Q - 



> 0. 



where; II = expected excess of market price of r, over future 
extraction costs of resource r in field z. 

= expected excess of market price of r, over future 



rz 



import costs of resource r in field z. 



n rz = expected excess of market price of r, over 
capital costs. 

Q r = expected demand for resource r. 



A set of equations paralleling (l)-(8) is thus required 
for "new fields" which result from exploration expenditures. 
The relations have the same interpretations as (l)-(8) but 
simply apply to field z=l,...,Z. 

(11) X t+1 = X t -W 1 +h t 

' rz rz rz rz 

(12) w* = <j> (n* , p t , Q 1 ) 

v ' rz wz rz rz r 



O') Q* = ^(Pj. P 2 P R ) 

( 13 ) 

(14) S t+1 = S t +W t +M t -A t 

v ' rz rz rz rz rz 

< 15 > 4z = 4=1, 

< 16 > “Jz “ +mz<^z=l,...Z> pJz-1, 

(17) S t < N* (K* ) 

rz — rz rz ' 



t 

. . Z ’ y rz=l , . 


. .Z’ 


Q t ) 
r 


t 

. .Z’ y rz=l, . 


. .Z’ 


q" ) 

r ' 



for 
all 
z, r 
and 
t . 
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Note that in (11), an initial stock X® is not given as in 

rz 

(1). If expenditures e t are made during the same period t, 

X* z 



equation (11) becomes: 

(ll 1 ) X t+1 = h t -W* and in t+1 
\ ' rz rz 



rz 



(ll 2 ) X t+2 = X t+ 1 -w t+ 1 +h t+1 
rz rz rz rz 

and a "new field" z is brought into the model. 

There are two general types of capital stocks which are 

used in the system. Those used for the extraction of resources 

in each field j and z, , and those stocks used for the 

z z 

management and storage of stockpiles of each r in each field 

j and z, . In each field j and z, K 1 ? and are lxN 

J ’ r j ’ rz ° ’ j z 

vectors, where K . , K L , n=l,2,...,N are intended to reflect, 

nj nz 

first a quantitative measure of various specific capital 
items (pumps, buildings, machinery, etc.) and second, a 
qualitative measure of capital items in terms of technologies. 
Thus, . may be a petroleum refinery with hydrocracking 

J 

capabilities, and a physical plant identical to 

except that hydrocracking capabilities do not exist. For any 

n, n, marginal investments, substantively different from 

normal capital-formation notion of investment (to be denoted 

as , later in this model), may have the effect of 

3 z 

"converting" . capital stocks to . capital stocks, where 
nj nj 

• 

such "conversions" have the effect of "converting" K_ . capi- 

II O 

tal stocks to K*. capital stocks, where such "conversion" 

nj 

have the effect of changing the technological structure of 
the extraction processing process. 
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Let o- . be the investment expenditures which converts 
n nj 

K . capital stocks to K_ . stocks, o ^ to convert K- . to 

nj nj n nj nj 

K . , o- the investment expenditures which converts K 
nj n nz nz 

capital stocks to K_ stocks, and a - to convert K- to K 

nz ’ n nz nz nz 

capital stocks. Assuming the conversion expenditures to have 
the functional form. 



( 18 ) 

( 19 ) 

( 20 ) 
( 21 ) 



O- . = <P -r(ru, tv , Q t ) 
n nj T aj v j’ 



o t - . = <j> .(n*, n t , Q t ) 
n nj T aj j’ 



J 

t 

'j 



0- = <j> ~(n t , n t , Q t ) 

n nz T az z z r 

0^ - = <b (n t , ) 

n nz T 0z v z z ^r' 



R 



where; II. = Z (IT.), n = Z (n ), 



r=l 



rz‘ 



J ~=i' r J 

= z (n^,) , rv = £ (tv ) • 



r=l 



r J 



r=l 



rz' 



If f=tj)j(0) is the quantitative measure of the resulting changes 
in the capital stocks after conversion expenditure 0 , then 



(22) 


f* 

n 


■ (b „-(0- 

nj T fj v n 


nj 


(23) 


f t 

n 


nj “ ^f j ^ °n 


n j 


(24) 


f 5 

n 


, t 

= <f> -(0- 
nz T fz n 


nz 


(25) 


f t 

n 


_ , ,t 
= <b„ (0 
nz r fz n 


nz 



Define V as additional investme. • i to capital stocks 
which are used for the extraction of resources in each field, 
and as additional investments to capital stocks which are 
used for the management, import and storage of stockpiles, 
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and having the general formats: 



(26) 


V*. = 
nj 


<j> .(n*, 

VJ V J 


4 


t t 

V V 


Q*) 

r 




(27) 


V* = 
nz 


<l> (n 1 , 

T VZ v z’ 


p z’ 


t t 

V V 


Q*) 

r 




(28) 


v t . = 

rj 


.(n t 

Vj v r j 


P rj 


t 

’ y rj ’ 


t 

n r j ’ 


«4> 


(29) 


V t = 
rz 


<S>- (n 1 , 
Y vz rz 


4 


’ y rz’ 


t 

n rz’ 


<£> 


Where 


t 

* P J = 


R t 

= E (p ■. 

r=l rJ 


), 









t _ R . t . 

y,- = £ (y ri ). 

J r=1 rj 



P ^ = A ( ^- 



and vt = S (y* ). 
z r=l rz 



Also define in general format: 

(30) D 1 . = .(W t ., V t ., K 1 . ) 

v ' nj Y dj v rj* nj nj' 

(31) D t = 4 . (W t , V 1 , K 1 ) 

nz T dz rz nz nz 

(32) D 1 . = 4> ,_ ( V 1 . , K 1 . ) 

v rj dj v rj’ rj 

(33) D t = 4>^,-( V 1 ' , K l ) 

rz dz rz rz 

where D as a net depreciation function for extraction capital 

gjj gjg gjg 

stocks, with _> 0, -gy £ 0, g^ >. 0, and D as a net depreci- 
ation function of capital stocks used for management of 

3D 3D 

stockpiles with — _< 0, — 0. Then the following transi- 

9 V 9K 

tion equations are assumed to describe the time path for 



capital stocks: 
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( 34 ) 



j— t-1 , . . . T . 

t=l, . . . ,T. 
t = l T 

v J- y m 9 m y y 

t=l T 

A word of explanation may be required at this point in 
order to clarify the meaning of (34)-(37) above. The first 
two of the transition equations specify that the capital 
stock of type n available for use in extraction field j 
(or z) at any point in time depends on the size of the stock 
in preceding period less net depreciation over the preceding 
period (which may actually be negative if investment in new 
stocks is large and/or if extractions are small), plus the 
net conversion of other types of capital stocks to type n. 
Thus, the latter summation term takes account of the amounts 
of other classes of extractive capital (n^n, 5=1,..., N) con- 
verted to type n, less that of type n, converted to other 
types. 

Equations (36) and (37) are required for the available 
capital used to manage and store the, up to R, different 
energy resources in field j (or z) at time t+1. Since it is 



K t+1 = K t .-D t .+ Z (f t ~.-f- .), 

nj nj nj fi=1 n nj n nj' 1 

n^n 

given, for all n, 5=1,..., N; 
N 

(35) K t + 1 = K 1 ' -D t + Z (f t _ -f- ), 

nz nz nz - n nz n nz 

n=l 

fi^n 

for all n, 5=1,..., N; j=l,...,J; 



(36) K t+1 = ,-D t . , for all j=l,...,J; 

v rj rj rj 

tfO 

K . given 
rj 

(37) K t+1 = K t -D t , for all z=l,...,Z; 

v 7 rz rz rz’ > > > 
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essentially assumed homogeneity of this class of capital for 
any given type of energy resource, there are no conversions 
to quantify. 

Next define as the quantity of energy resource r 
available for use in the production of "final energy commod- 
ities" at time t. 

Thus uVmay be thought of as the volume of petroleum 

(crude) made available to refineries, coal at trackside, etc. 

The term U t then includes the sum of draw downs from 
r 

stockpiles of r from all resource fields, j, and z, for cur- 
rent consumption, then; 

t J t Z t 

(38) U = Z A .+ E A c , for all r=l,...,R; t=l,...,T. 
r j=l rJ z=l rz 

Next it is assumed, there exist L different "final energy 

"t t h 

commodities" and define as the amount of the 1 type 
produced at time t, (1=1,..., L). The amount of each produced 
depends upon the amount of the R possible energy resources 
used as inputs "U^,", as well as up to Q different kinds of 

A f 

capital stocks "K used in the production process. Thus 
the production functions could be written as: 



(39) 



with 



A i = v u £i> • • • > u 



Rl ; K ll’ 



’Kqi>’ 



for all 1=1, 



,L; 



9A* 

— — ^ 0; -~zr >_ 0 for all r, q. 
9U^ 9iq 



t = l T 

1/ X y • • • y X • 



Note that obviously all energy resources might not be used as 

9 A 

inputs in each production function process; then ^ — may be 

° U rl 

zero for some r in each 1. 
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Now let each unit of energy commodity 1 contain BTU's. 
If E* equals the minimum amount of BTU required in period t, 
then the problem below must be constrained by: 

L t t 

(40) Z ,A^ > E , for all t=l,...,T. 

1=1 1 1 

Let >. 0 equal the minimum required quantity of the 1 L 
specific energy commodity in period t. This implies. 

(41) A* > H*, for all 1=1, . . . ,L; t=l, . . . ,T. 

C. THE CRITERION FUNCTION 

The energy independence policy presented by President 
Nixon has been used as an application of the model introduced 
in previous section and for the formulation of its criterion 
function . 

The basic concept of this policy arises from the threats 
to the security of the nation because of foreign supply 
interruptions, and the assumption that the threats most likely 
to be encountered, are best met through domestic self-suffi- 
ciency. Therefore, it calls for more storage of resources, 
their higher production, and expansion of explorations for 
new deposits and sources of energy, aimed toward the realiza- 
tion of this concept by period T. (1980) Within the decision 
situation described by the conditions and restrictions (1)- 
(41) of preceding section, concern here is with the indirect 
selection of values for extraction rates of natural resources, 
the expenditure rate for exploration, rate of investment of 
capital stocks, rate of production of energy commodities to 
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respond to energy requirements and production such that the 
following expression, the present value of total capital out- 
lays and consumer costs over "T" periods of time be minimized 
T 

(42) Min E C t 3 t +<K r (X T ,K T ) 3 t 
t=l 1 

Subject to (1)— (41) 

where; C*' = Sum of all capital outlays and consumer costs. 

= A concave terminal value function for all the 
natural resources and capital stocks. [23] 

0 = l/(l+r), where r is the appropriate positive 

discount rate. Reflects the discounting factor. 
The objective function in a detailed form could be written as 
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t=l 
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+E o- 
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nz 
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n nj 
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+E. 


wt . 


•C Wt 0t+E 


w t . 


c wt 


gt + E. 
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rj 
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rj rz 


rz 
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rj 


rj rj 




+ E 
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• C S t gt + E. 


wt 

M . • 


C mt 6 


t + E M 
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•C mt 6 t 
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rz 


rz 


rz rj 


rj 


rj 


rz 


rz 


rz 




,(x T 


T T 
,K i )0 i . 















Subject to equations (1)-(41) 



R Z R J 

where; E = E E ; E. = E E; E 



rz 



r=l z=l 



rj 



r=l j=l nj 



T T T 

X = E. X . + Z X 1 
rj rj rz rz 



np rp rp m rp a rp 

K = EK +EK +E. K 1 . + E K +E K n . 
j j z z r j rj rz rz i q lq 



N J 
_E E 
n=l j=l 
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D. EQUILIBRIUM CONDITIONS 

The optimization problem expressed by (43) and equations 
(1) through (41) could be solved by Lagrangian method. As- 
suming that the functions in (1) through (43) are continuous 
and possess first partial derivatives with respect to each 
argument, the following Lagrangian expression may be formed. 



(44) L = E (E e l g t +E. V t ,g t +Z V 1 B t +S. V t .B t +S V 1 
v rz rz nj nj nz nz rj rj rz rz 

+ E a t _.B t +E a t _ B t +2. a- .& t +Z a- 
_ . n nj - • n nz nj n nj nz n nz 

11 J 11 J 

+ E. W t .•C W ^B t +S W* •C wt B t +S. S* . •C st e t +E S t •C St B t 
rj rj rj rz rz rz rj rj rj rz rz rz 

+ E. M t . •C m ^B t +S •c mt B t - 1 x t • [X t .-X° .+ E W t .] 
rj rj rj rz rz rz r j 1 rj rj ^ r j 



- E. 2 X ^ [\T .-4> ^ . ( n 



rj 



t 

^wj 



t 

rj 



rj rj 

p 2 , 



- E. 4 X t 



rj 



. [W* .~g t .( X t . 
r j rj rj rj 



pJj, <4)i 



••Vi 



Kj), 



- 5 X t . [S t .-S° E W* E M t .+ E A t . ] 
r J rj r j T=1 rj t=1 rj T=1 r j 

-e. 6 x t .[A t .-4>t .(n t . , T ,P t ■ 
rj r j rj r 6j v r j=l , . . . J ’ K r j=l , . . 



J’ P rj=l, 






- E. 7 X t . [M* .-(p 1 . (n t . 1 T ,p t . n T ,y t . , 
rj r j rj y mj rj=l , . . J H r j=l , . . J H r j=l , . . J r 

] 



- E. 8 X t . [S t .-N* .(K l - ) ] - 9 A 1 



rj rj 



E 10 X t 



,t 

*rj 

t 



[h t (e t )-X 
rz rz rz rz 



rz 



[e~ -<p (n t ,p t .n* ,Q t ) ] •B t 

rz rz e rz rz rz r 



- 1 1 X t [X t + E W 1 - E h t ] 
rz rz rz rz J 

T =1 T =1 



. y l 2 x t 

rz A rz 



A v [W t -<p (n* 

rz r wz v 



rz'^z' Q r>> 



-E 13 X t [W t - g t (X t , K*)] 
rz rz 1 rz °rz rz z 1 
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_ 14 x t [s t _ w t _ M t + £ A t ] 

rz rz , r _ 1 rz i rz 1 rz J 
T-l T =1 T =1 

* ‘ 5X L |4 L-*L (n L=l, . .Z' p rz=l, . .z-^z-1, . .Z’ Q r>l 



iFz X rz [M rz <J) mz (n rz=l, . .Z ,p rz=l, . .Z’ u rz=l, . .Z ,Q r ) ] 
-JL I7 x!L[s! -N* (K* )]-E, le X^.B t [a^ .-4 ,(n^,n^,Q!) ] 



rz rz 1 rz rz rz 7 nj nj’ 
19 X^.fa t _.-<*> .rn^.Tij ,q£) ]B t 



n nj T aj 



-E 

nj 



-E 

nz 



2 0 ■ 



nz 



nj 



nj 
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nj 


oj J 
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r t 

[a- 




■4> .(n* 


nz 


L n 


nz 


r o 2 V z 
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r t 

[a 




■4> (n t 


nz 


n 


nz 


r az z 


t 

nj 


■*s 


nj" 








nj 



l n j 1 n nj v fj v n nj 



-E 24 X t [f^ -4>„_(a^ ) ] -E 25 X- [f t - -<()- (a t - )] 

nz nz n nz r fz n nz - nz 1 n nz fz n nz J 

nz 

-e. 26 x t .e t [v t .- 4 ) .(n t ,p t ,u t ,n^,Q t )] 

nj nj M nj Y vj v j j j j r /J 



-x”x^eXz-* V z< n z-P^V ll z'‘ 3 z>] 



nz 



nz 



-E. 28 x t .B t [v t .-<j)_ .(n t . , p t . ,vi t . ,n t ■ ,Q t ) ] 
rj rj 1 rj y vj v r j ’ K r j ’ M r j rj’^r /J 

-E 6 [V -d>- (n ,p ,y , n ,Q )] 

rz rz 1 rz y vz v rz ,K rz’ rz* rz r J 

-E. 30 X t . [D* .-4> . .(W* . , V t ., K* . > J 



nj n j nj Y dj v r j ' nj ' nj 

- ;l>, *L id L-W' ? L- v -- K ^~>! 



nz 



nj 

nz 



nz 



-E. 32 X t . [D t .-4) ,_(V t . ,K t .) ] 
rj r J 1 rj Y dj v r j r j 

-E 33 X t [D 1 -4)._(V t , K 1 )] 
rz rz rz dz rz rz J 



t f; N 

_ 3 4 ■> t rT ,t T ,0 ; ” /Js t 



nj 1 nj nj T=1 nj T=1 - =1 ' n nj n nj 



- 35 X t 



x;,[k:,-k„,- e d;.- e e (f; s ,-f£ r ,)) 

)] 



t ^ t ^ ^ t t 

[K l - E D l - E E (f l - -f_ 

nz nz^-nz^-.-^nnz n nz 
x=l t= 1 n=l 
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- 36 A t . [K l ,-K° .+ Z D t .]- 37 A t [K^ + Z D 1 ] } 



r J rj rj T=1 rj 



- rz rz _ - rz 

T = 1 



T 

Z 

t=T 



+ Z_{- 38 A^[A^-4) 1 (U^ 1 ,...U^ 1 ; , . . . Kq X ) ] 



- 39 A^[^Z^a 1 A^-E t ]- 40 A^[A^-Hj] }+ ^ T (X T ,K T )B t . 

The necessary conditions for minimizing equation (44) 
include the following equations, which are then respectively 
interpreted for their economic implications on the objective 
function. 

a<J t 

(45) = 3 t - 2 6 A^ 3 0 A^ . ( — 

r l n j n j gy l 



3V' 



nj 



nj 



4 * 

In equation (45) the term 26 A\g represents the discounted 
value of the incremental change in objective function due to 
a unit of change in the investment of capital stocks used 
for extraction of resource r in field j, defined by constrain- 

t 3<f) di 

ing equation (26). The terms 30 ^ n -[ — ^ ] represents the 



3 V 



nj 



implicit value of a marginal change in net depreciation of 
these capital stocks due to a marginal change in relative in- 
vestments. The necessary condition for optimality requires 



that 



3L 

av 1 . 

nj 



= 0, or 



. 3 4> j . . . . 

(45-1) 3 0 A^.[ % 1 ] = P - 2 6 A*| .p . ' 

nj 

A similar interpretation can be made for equations (46)- 
(52) below. 
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(46) 



9L 



0 



(47) 



(48) 



(49) 



(50) 



( 51 ) 



3 d) ^ 

— — = e t _ 27 x t 6 t - 3 1 A t [ — = 

av t nz ■ nz av* 

nz nz 



9<t> 



1 ^ "t r ^ _ q t 2 7 t p"t 

nz r 



(46-1) 3l A" [ — = 6"- zy X" 6 



nz 3v t 
nz 



3L 

av* . 

rj 



3 (J) 

= b 1 - 2 8 x t . g 1 - 3 2 x t . [ — Ji] = o 

rj rj^t 

rj 



32 T^ f 1 = 



3 +dj 

r*. 

rj 



(47-1) Jz x; A-zf-] = 3 - ZB x; b 

i j av 1 j 



Hi: 



4 - = B *-' 2 * 4 [— 4 ) = 0 

— t rz rz 



3V 



rz 



3V’ 



rz 



, a 4> j~ . , . 

(48-1) 3 3 X^ z [-^|5] = B t - 2 9 X^ z B t 

3 V 

rz 



3L 



3a* - . 
n nj 



, . , , a <f>3 • 

= B 1 - 1 9 xl ^ 6 t + 2 3 xl ^ [ — p — ] 



nj 



3a t -. 
n nj 



, 3(f) ~ • , , , 

23^t r fj j _ l 9 X t .g t -B t 



(49-!) ZJ x£ [ 

J 3a 



nj 



n nj 



3L 



aa t - 
n nz 



= 0 



. 3<f>^ 

25 x L ( d^> 

3a - 
n nz 



B t [ 2 1 X^ -1] 
nz 1 



3L 



3a^ 

n nj 



= 0 



f 3 4> • 

(Si-i) 22 X 4 (-4 
nj 3ai 



] = 






n nj 



0 
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(52) 



0 



3L = 

3ai 

n nz t 

(52-1) 



So- 
il nz 






(53) 



3L 



34>* . 



= 6 t *C W ^- 1 X t .- 2 X t .-“X t . + 5 X t .- 3 °X t . [ — = 0 
rj rj rj rj rj nj 1 aw t 



3W . 
rj 



*t W"t 

In equation (53) the term g *C . represents the discounted 

J 

t 9< ^di 

value of the extraction cost, and 30 X .[ — -^1 represents the 

"J SwV 

rj 

implicit value of the marginal change in net depreciation of 
the capital stocks due to a marginal change in extraction 
rate of resource r, where J X^j represents the value of the 
incremental change in objective function due to a unit of 
change in the extraction rate restricted by equation (1) 
and 2 X t . for equation (2) and 4 X t . for equation (4) and 5 X t . 

J -*■ J J 

for equation (5). The necessary condition for optimality 



requires that 



34> 



(53-1) 3 0 X t .[ — = g t *C W ^-( 1 X t . + 2 X t . + 4 X t .- 5 X t . ) 



nj^ w t 
rj 



rj rj rj rj rj 



A similar interpretation could be made for the equation 
(54) below. 

34> 



(54) 



3L 



3W 



= g t -c wt - 11 x t - 12 X t - 13 X t + 1 “X t - 31 X t [ — §2] = 0 



rz 



rz rz rz rz 



rz 



rz ‘3W l 
rz 



3 (j> 

(54-1) 31 X t [ — £2] = g t -C Wt -( 11 X t + 12 X t t 13 ! 1 - 14 X t ) 
v ' rz‘ .„t J rz v rz rz rz rz' 

rz 

9h t 

(55) = g^X* [ — B t + 11 x^_ = 0 



3e 
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rz ‘ se L 

rz 



rz 



rz 
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. 31i 

q X> rz 

In equation (55) the term 9 X [ — — ] represents the implicit 

rZ 3e* 
rz 

value of a marginal change in units discovered of the re- 
source r in field z due to a marginal change in the explora- 
tion expenditure restricted by equation (9). The term 
10 represents the discounted value of the incremental 
change in objective function due to a unit change in explora- 
tion expenditure described by equation (10). Similarly the 
term 11 represents the value of an incremental change in 
objective function due to a unit change in the quantity dis- 
covered of resource r described by equation (11). The 
necessary condition for optimality requires that 



. 3h 

( 55-1 ) 9 ^ z [ t 



3e 



rz ] = g t (i- io A; z )+ 



1 0 1 "t >4.11 



rz 



rz 



The most important consideration in this study, is the impact 
of price changes, on the objective function and other vari- 
ables in the model. This is described by (56), which is the 
partial derivative of equation (44) with respect to P^, the 
price of resource r in period t, 
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written 



,2 9->t p t 3 **rz . 3<t> vz 3p rz ^vz 3p rz 
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In equation (56-1) the terms outside of the brackets repre- 
sent the marginal change in the excess of markets price of 
resource r over different costs in different fields due to 
a marginal change in the market price of resource r. And 
the terms inside the brackets (some of them discounted) 
represent the implicit value of the marginal changes, in 
different variables, influenced by the price of resource r, 
due to a marginal change in the excess price over the costs. 

In optimality the total sum of these terms in the left 
hand side of equation (56-1) must be equal to the term on 



* 3<f> 

~ t~ Y q * 

the right hand side - A [ — J] which represents the implicit 

r 



value of the marginal change in units demanded of resource 



54 



r, due to a marginal change in its market price, the demand 
function is described by equation (3). 

The above mentioned conditions simultaneously solved 
with the following equations would determine the values of 
multipliers and the optimum rates of exploration, investment, 
extraction, importation and conversion expenditure. 



(57) 





(58) -^r 




rz 




rz 



0 



(58-1) 16 A* - 



rz 







(60-1) 1 4 A 1 + 1 7 A t 



rz 




(61) = 0; for all As. 
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V. CONCLUSIONS 



An attempt has been made to bring together and formalize 
some of the factors that seem to exist in "Supply and Demand 
for Energy". After a brief review of the world energy con- 
sumption it is shown that there is a broad parallel between 
economic growth and energy consumption. It is important to 
explore the relationship as precisely as possible and thus 
more intensive research is needed to determine: 

a. The impact on energy consumption of technological 
efficiencies in consuming countries. 

b. The "average" nation-wide levels of energy consump- 
tion required for given level of GNP. 

c. The effect of above mentioned topics studied in terms 
of specific energy sources. In addition, the substi- 
tutability of one energy source for another in satis- 
fying requirement should be considered. 

d. A more specific effect of relative prices of resources 
on energy consumption. 

e. Impacts of environmental controls and pollution 
abatement upon the volume of energy demanded and upon 
shifts to substitute energy sources. 

f. Impacts of economic and financial policies such as 
balance of payments or balance of trade considera- 
tions upon the volume of energy consumption. 

By introducing the above mentioned considerations in the 
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empirical analysis of past data, a greater insight might be 
gained into future projection of energy requirements. 

In the analytical section of this study the Cummings and 
Whipple model was extended by posting a relationship between 
extraction, storage, import, exploration and market prices 
of resources and relative cost factors. 

The model does not include social values and environmen- 
tal considerations. 

The model would allow for price regulatory policy deci- 
sions, and considers both supplies and demands of resources, 
and their interactions is achieved through the dynamic mar- 
ket mechanism. It is designed for repeated use and consid- 
ers the technological impact and conversion of capital stocks 

The model is disaggregated by regions, covers all energy 
forms and takes into account the interfuel substitution 
brought about by price changes. 

It is believed that expansion of the model to cover the 
following additional areas, would make it more useful in con- 
junctive management. 

i 

1) Inclusion of environmental considerations and pollu- 
tion abatement. 

2) Inclusion of specific functional forms for extraction 
import, investment, etc. These could be obtained by statis- 
tic estimation. 

3) The inclusion of uncertainty in exploration expendi- 
tures. 

4) Inclusion of transportation problems. 
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5) Refinement of the model so that a computer simulation 



might be employed. 
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APPENDIX A 



1 . Compound interest and rate of growth 

If $y 0 is invested at compound interest at r% per year, 
compounded yearly then the amount after t years is $y where: 

(1) y = yo(l + r) t 

if the interest is added n times a year, then 

r nt 

(2) y = y o(l + 

Discontinuity is an essential feature of this compound inter 
est problem. The more frequently is interest added, the 
larger is the amount of a given sum at the end of any period 
Letting n, the number of the frequency of compounding 
interest to take larger and larger values, it is clear that 

1 n 

(1 + — ) tends to a definite limit and it can be shown that 

1 n 

e = lim (1 + -) ; e = 2.71828 
n 

n-*-eo 

Returning to the case of compound interest 

nt n_ rt 1 rt 

y = yo(l + J) = yo((l + ^)r} = y 0 ((l + ^) m > 

1 m 

for n-*- oo : (1 + — ) = e then 

’ v m 

(3) y = y 0 e 

We have now derived a concept for continuous interest. 

Note: This model represents a constant percent rate of 

growth (interest rate). It is clear that the percent rate 
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of growth at time t is r(t) = ■ y - ( 100 ) where ^ represents 
the rate of change. 



2. Present Values 

a. If the interest is added once yearly at r% then y 
is the present value of y available after t years; from Eq. 
( 1 ): 

(4) y o = — - — r 

( 1+rT 

b. If interest is added n times a year at r% from Eq. 

( 2 ): 

<5 > y - — frt 

o+f) 

c. Finally, if interest is added continuously at r% 
from Eq. (3) it implies 



(6) y o = -Jjr or y 0 = y*e rt 
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APPENDIX B 



Calculus of variations - Euler's Equation . 
If the extreme values of the integral 



1 



u = / f (t ,q(t) ,^)dt 



are required for all possible variation in the function 
q = q(t) such that q(t 0 )=qo and q(t 1 )=q J where (t 0 ,q 0 ) and 
(tj,qj) are fixed points. The function 



f(t) E f(t,q,^) E f{t,q(t),q(t)} 



which gives variable u on integration, depends on the variabl 
t, on the variable function q=q(t) and on the first deriva- 
tive q(t). This is a function of t given in the function of 
functions form. In solving the problems, the function q(t) 
is taking in the form q=q(t ; a, 3, y, . . . ) where q(t ; a, 3, y, . . . ) 
is assumed to be a fixed function (with a continuous deriva- 
tive) and where a,3,y,... are parameters. Given differential 
increments 6a, 63, 6y, ... to the parameters, we derive the 

corresponding variations 6q and 6q in the function q and its 

derivative q = ^ : 



(1) 

( 2 ) 



6q 

sq 



|a Sa + |a + |a 6y + 

8a 83 8y 1 



and 



- «H> - &&>*« + + 

= + ft ( le )5B + + ••• 



dt v 33' up ' dt v 3y 

8ci „ ^ d 

w 1 ' 



dt v 8a 

= ft ( l^ 6a + 83 66 + l? 6Y + ‘ ‘ ° " '3t (6q) 



1 



Reference 7, Page 521. 



61 



The symbol "6" is used in order to distinguish parametric 
differentials from "d" deferring to variation in variable t. 
The function f(t,q,^-) and the integral u can now be consid- 
ered as dependent on the parameters a,B,y, ... and the varia- 
tion in their values are obtained as: 



_ 3f,„ . 9f x • 

6f 3q 6Q 9q 6q 

from Eq. (2) results 

/on x-p 3f x , 3f d , x v 

(S) « = + 14 dt ({q) 



and 

6u = 

From Eq. (3) 



1 1 

/ f ( t , q , q )dt 
to 



1 1 

/ 6f dt 
to 



(4) 



du - /(jfdq)dt + / { H^(d q)) dt 

t 0 



. V ,8f d 

to i- 0 

and finally be appropriate substitution and rearrangement. 
The expression for variation in u becomes: 

. ,3f. , t *. V ,3f d ,3f., r .. 

di = t^dq) + / {jj - dt<34)}« q dt. 

L 0 X 0 

Since by boundary conditions, the curve q=q(t) always passes 
through two fixed points t=t 0 and t=t i , it follows that the 
6q=0 at these points. Hence 



of 1 1 

[fgdqj = 0 

t 0 



and so 



(5) du = /‘(f| - |f(|f)}dq dt 

t 0 
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The necessary condition for the u to have an extreme value 
for variation in the function q=q(t) as obtained by varying 
the parameters a,B,y,... is that Su=0 for all values of 6a, 
6B, 6y, i.e. for all values of 6q from Eq. (5). This 

is only true if 



(6) 9q dt ( 3q ) " 0 



This result, known as Euler's equation, determines that 
function q=q(t) which maximizes or minimizes the value of u. 
It is a differential equation which must be solved to give 
the function q=q(t) sought. Euler's equation, however, is 
only a necessary condition for extreme values of u. For 
distinguishing between maximum and minimum values, a general 
criterion is not readily obtainable, in simple practical 
case, it can usually be judged from general reasoning. 
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APPENDIX C 



OECD Countries used in the Analysis 

1. Canada 

2. United States 

3. Japan 

4. Australia *** 

5. Austria 

6. Belgium 

7. Denmark * 

8. Finland * 

9. France 

10. Germany 

11. Greece * 

12. Iceland ** * Not used in Natural Gas 

Analysis . 

13. Ireland * 

** Used only in Aviation Fuel, 

14. Italy Motor Gasoline, Kerosine 

and Gas Oil Analysis. 

15. Netherlands 

*** Used only in Energy /Capita 

16. Norway * Analysis except for the 

years 1966, 1967. 

17. Portugal * 

18. Spain 

19 . Sweden * 

20. Switzerland * 

21. Turkey * 

22. United Kingdom 
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ENERGY/ CAPITA AND GNP/ CAPITA 
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Figure 2. ENERGY/CAPITA AND GNP/CAPITA 
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Figure 3. CRUDE PETROLEUM/ CAPITA AND GNP/CAPITA 
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Figure 4. CRUDE PETROLEUM/CAPITA AND GNP/CAPITA 
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Figure 5. MOTOR GASOLINE/CAPITA AND GNP/CAPITA 
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Figure 6. MOTOR GASOLINE /CAPITA AND GNP/CAPITA 
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Figure 7. AVIATION FUEL/CAPITA AND GNP/CAPITA 
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Figure 8. AVIATION FUEL/CAPITA AND GNP/CAPITA 
Constant Elasticity Model 
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Figure 9. HARD COAL/CAPITA AND GNP/CAPITA 
Linear Model 
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Figure 10. HARD COAL/CAPITA AND GNP/CAPITA 
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Figure 11. KEROS I NE/ CAPITA AND GNP/CAPITA 
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Figure 12. KEROS I NE/ CAPITA AND GNP/CAPITA 
Constant Elasticity Model 
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13. GAS OIL/CAPITA AND GNP/CAPITA 
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Figure 14. GAS OIL/CAPITA AND GNP/CAPITA 
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Figure 15. NATURAL GAS/CAPITA 

Linear Model 
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Figure 16. NATURAL GAS/CAP fTA AND GNP/CAPITA 
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